How to Define Variation of Physical Properties Normal to an Undulating 

One-Dimensional Object 
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One-dimensional flexible objects are abundant in physics, from polymers to vortex lines to defect 
lines and many more. These objects structure their environment and it is natural to assume that 
the influence these objects exert on their environment depends on the distance from the line-object. 
But how should this be defined? We argue here that there is an intrinsic length scale along the 
undulating line that is a measure of its "stiffness" (i.e., orientational persistence), which yields a 
natural way of defining the variation of physical properties normal to the undulating line. We 
exemplify how this normal variation can be determined from a computer simulation for the case of 
a so-called bottle-brush polymer, where side chains are grafted onto a flexible backbone. 

PACS numbers: 82.35.Lr,02.70.Uu 



One dimensional objects are characteristic for the 
physics of quite diverse phenomena. Most prominent 
may be linearpolymers [3, Hj] , be they synthetic ones or 
biopolymers [3j, and vortex matter in high T c supercon- 
ductors [1,01 • But line objects also determine the char- 
acteristics of the structure of (poly-)crystals in the form 
of defect lines [6( and they were argued to determine the 
structure of metallic glasses [7j in the form of disclination 
lines. This list, of course, makes no claim to complete- 
ness. In all cases, however, physical properties around 
the line objects will change depending on distance to the 
object. For polyelectrolyte chains like, e.g., DNA, one is 
for instance interested in the distribution of counterions 
around the chain [f| and the phenomenon of counterion 
condensation on the (linear) chain. In the case of vortex 
lines, the order parameter distribution around these lines 
is important Q , and there exists also a rich phenomenol- 
ogy of interactions between vortex matter and the defect 
structure of the underlying lattice [lfjj. It is also well 
known that dislocation lines in metals control their plas- 
tic deformation. Defect lines in ordered materials and the 
strain fields around them are directly observable in col- 
loidal crystals and nematic materials ll], 12, IH and are 
of importance for instance for the performance of pho- 
tonic crystals [1J|. In all these cases, physical quantities 
can be expected to change in space according to the dis- 
tance of a given point from the line object. Such a normal 
or radial distance is clearly defined for a straight line, and 
sometimes one uses models assuming infinite rigidity to 
look at radial variation, like for example the cell model 
of polyelectrolytes , or chain molecules tethered to a 
straight line as a model for a molecular bottlebrush , 
etc. However, the usefulness of all these models has been 
completely uncertain. In reality, however, these line ob- 
jects are not infinitely rigid and thus they undulate due 
to thermal fluctuations. This undulation leads to an in- 
teraction between different parts of the lines through the 
structuring effect they have on their environment. The 
question of how to define a radial distance from the line 



object in this case naturally arises when one studies such 
physical systems with computer simulations, but also for 
experiments [l7j |. when the structure of a system can be 
observed for instance microscopically or tomographically. 

From the differential geometry of space curves we have 
an easy answer to this problem locally. Just calculate 
the two normal vectors (in 3d) to the line. And globally, 
when the line object is curled up to a globule or coil, 
there is no natural normal direction. While the latter is 
true also when we take into account the granularity of 
matter, the local definition seizes to be helpful on the 
atomic scale. Here we will find preferred angles between 
consecutive segments of the line given either by the un- 
derlying chemical nature of the object like in polymers or 
by an underlying lattice structure like in vortex matter or 
defect lines. The natural length scale for a definition of 
normal variation will therefore lie in between this local - 
and very system specific - and the global - in many cases 
isotropic - scale. We will discuss in the following how an 
optimal choice for this intermediate length scale can be 
identified both in simulations and in experiments provid- 
ing real space structural information. The procedure will 
be completely general, however, we will exemplify it for 
the case of a so-called bottle brush polymer [lq . 

Our line object will therefore locally be made out of 
straight segments - we call them bonds from now on - 
making up the linear backbone of a bottle brush polymer 
shown in Fig. Q] We regard a segment of the backbone 
of this bottle brush made out of l c consecutive bonds. 
For this segment we can define an equivalent cylinder 
axis as follows: each unit bond vector can be seen 
as the normal to a plane, in which we define cylinder 
coordinates. The average orientation of these planes is 
given as 
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This definition of n is similar to how one would define the 
normal to an undulating surface. When we assume a tri- 
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defining what is meant by variation normal to the back- 
bone of the bottle-brush polymer. Our construction en- 
sures that each monomer is contributing not more than 
once to the determination of the radial density profile; 
however near the backbone chain ends some monomers 
may not be counted at all (when their z-coordinates do 
not fall into the allowed range for any local cylindrical 
coordinate system). This is an example of chain end ef- 
fects occurring whenever the undulating line has only a 
finite length. The normal variation defined in this way 



FIG. 1: Snapshot of a bottle brush polymer simulated us- 
ing the bond fluctuation model. The backbone consists of 
Nt — 387 repeat units and onto each of these a side chain of 
length N = 48 is grafted. One segment of the backbone con- 
sisting of l c bonds is plotted next to the chain. Such segments 
define a local coordinate system also indicated in the figure 
and cylinders surrounding the simulated chain as shown (see 
text). 



angulation of the surface, each triangle defines a unique 
normal vector, and the average over these normal vec- 
tors for adjacent triangles defines the surface normal on a 
coarse-grained scale. Using h we can define an equivalent 
straight cylinder for our undulating line on the coarse- 
grained scale l c . The z-axis is given by h, it goes through 
the center of mass (cm.) of the backbone segment made 
of l c + 1 monomers and has a length given by the end-to- 
end distance R% = \\r n+ i c — f n \\ of the backbone segment 
between monomers n and n + l c . In the normal plane 
to this axis we employ cylinder coordinates to define the 
radial distance of a given point to the axis of the cylin- 
der. Each point in space is considered to count for the 
normal variation within this cylinder segment, when its 
z-coordinate in the local coordinate frame lies between 
z = and z = R c e . By construction, each point in space 
may lie in more than one of these cylindrical segments, 
in which case it is, however, only counted at the small- 
est distance to these segments (when degeneracy occurs 
where a monomer has the same distance to the axis in 
two cylinder segments, it should be counted for both seg- 
ments with weight 1/2). Evaluating whichever scalar or 
tensorial physical property is characteristic for the con- 
sidered system at these points, one can in this way map 
out a variation of this property normal to the undulating 
line for which one performed this construction. 

In our case of bottle brush polymers we study the 
normal variation of the density of side chain monomers 
with respect to the backbone of the polymer. This den- 
sity variation is something one t ypi call y tries to extract 
from scattering experiments ID, 1^, [2(J, Hl[ . Due to un- 
avoidable approximations in the analysis of the scattering 
data, there has been, however, a controversy in the lit- 
he form of this spatial variation. 
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To resolve this controversy we apply our procedure of 



FIG. 2: Variation of the density of side chain monomers nor- 
mal to the backbone of the bottle brush for Nt = 131 and 
N = 48. Results for several choices of the coarse-graining 
length lc along the backbone are shown. All lengths are mea- 
sured in units of the lattice spacing of the underlying simple 
cubic lattice. 



depends on the coarse-graining length scale l c along the 
undulating line. Figure [2] shows the dependence of the 
normal density on this length scale for a backbone of 
length Nf, = 131 and side chain length N = 48. The data 
are for a Monte Carlo simulation of the bond-fluctuation 
lattice model [22j]. Clearly, the form of the curves for 
very small (l c = 3) and very large (l c = 130) values of 
the coarse-graining length significantly differs from the 
behavior for intermediate scales. The behavior for small 
l c gives comparable results to an approach taken in [23| 
to measure the counterion concentration around a poly- 
electrolyte chain. We now want to argue that at the in- 
termediate scales there exists an (almost sharply defined) 
optimal length scale for coarse-graining. 

This length scale has to be connected with the intrin- 
sic stiffness of the undulating line object. In polymer 
physics this is called the persistence length of the poly- 
mer chain (here the backbone of the bottle brush). The 
persistence length is known to be the natural scale for a 
coarse-graining of intramolecular behavior, i.e., variation 
as a function of chemical distance along the chain. We 
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FIG. 3: Projection of the end-to-end vector of the back- 
bone of the bottle brush onto the local bond orientation as 
a function of position of the bond along the backbone. The 
backbone length is Nb = 259 and the side chain lengths are 
N — 6, 12, 24, 36 and 48 from bottom to top. Dotted lines 
indicate values of the persistence length. 



suggest here, that it is also a natural scale for construc- 
tion of a coarse grained description of property variation 
as a function of spatial distance to the chain. There 
are several prescriptions of how to obtain the persistence 
length from an ensemble of configurations of a chain. One 
definition uses the projection of the end-to-end vector of 
the chain onto the unit bond vector of bond k, measured 
in terms of the average length of this bond 

l p (k) = (b k ■ R e /\b k \ 2 ) . (2) 

In a plot of l p (k) as a function of position k along the 

chain as shown in Fig.[3]one can clearly identify a plateau 

regime. The height of this plateau defines the persistence 

length l^p \ Another definition employs the bond vector 

orientational correlation function (cos#(s)} — (e$ • ei+ s ) 

where the angular brackets include an average over i. For 

idealized polymer models this function decays exponen- 
ts 2) 

tially as a G exp{— s/l p } (the amplitude prefactor is nec- 
essary for discrete models for which the exponential does 
not extrapolate to one for s — > 0). In reality, this func- 
tion generally is not a single exponential, and we propose 
to define the persistence length for arbitrary undulating 
lines as 

, N b -1 

V = sjispnj £ < ™ 9<s) > ■ <3) 

v 7 s—l 

(2) 

If l p >> 1 and the decay of < cos 9(s) > with s is a sin- 
gle exponential, one can for Nb — > co transform the sum 
to an integral and obtains again the standard definition 



of Ij, . It must be noted, however, that for real poly- 
mer chains (which have excluded volume interactions) 
one finds power law decays rather than exponential de- 
cay for Nb — > oo[13, US Ell . For isolated self-avoiding 
walks (good solvent conditions) < cos 9(s) > oc s 2v ~ 2 
and then l p 2 ^ oc Zp 1 ' oc N 2 "' 1 2A\. In dense melts 

< cos 9(s) > oc s" 3 / 2 HH, l p 2) {Eq. (3)} then is domi- 
nated by the behavior for small s, and does not depend 
on Nb for lager Nb- Thus the use of Eq.(3) requires care. 
Both estimates of the persistence length agree within 10% 
with each other and show an increase by a factor of about 
4 (see Fig. [3]) when we increase the side chain length from 
N = 6 to N = 48. This is in excellent agreement with 
experiment [271 ] and agrees with simulations in (28j but is 
a stronger variation than found in [29(. For the present 
purposes the important point is that the stiffness length 
scale l p (and hence the coarse-graining length scale l c in 
Fig. [l} can be varied over a wide range, which allows us 
to assess the validity of our general concepts. 
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FIG. 4: Density variation of side chain monomers normal to 
an undulating backbone of length Nb = 131 for the choice 
l c ~ l P and for the rigid backbone. The four choices of side 
chain length are N = 6, 12, 24 and 48 from left to right. 



When we now take the persistence length as the nat- 
ural scale for the definition of the coarse-graining length 
in Fig. [1] i.e. l c = l p , we find that the normal vari- 
ation determined in this way is in excellent agreement 
with what one would determine for a completely rigid 
backbone. In Fig. 3] we compare the radial density vari- 
ation for a backbone length Nb = 131 and four differ- 
ent side chain lengths N — 6,12,24,48 obtained using 
l c = lp to the radial density variation around a rigid 
backbone [32j]. Obviously, our definition of normal vari- 
ation with respect to the undulating backbone for this 
choice of coarse-graining length nicely agrees with the 
variation around the rigid backbone at small to inter- 
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mediate distances. In this way, we justify the often used 
rigid backbone model for densely grafted bottle-brushes. 
For large distances, the density measured for a flexible 
backbone has to be larger than the one measured for the 
rigid backbone, as the former contains contributions from 
remote segments of the backbone which can bend back 
on itself. From our experience, we can vary l c by about 
20 % around the choice l c = l p without much affecting 
the agreement observed in Fig. 21 Therefore, the per- 
sistence length provides the natural length scale for the 
coarse-graining depicted in Fig. [T] Reinterpreting Fig. Q] 
from a different perspective, we can say that an analy- 
sis of the variation of physical properties normal to an 
infinitely stiff line-object gives a good approximation to 
the behavior for real, thermally fluctuating lines on the 
scale of the persistence length. 

To conclude, we have outlined a constructive procedure 
how to define and measure normal variation of physical 
properties with respect to an undulating line. This prob- 
lem is relevant for a broad range of physical situations 
ranging from linear polymers, to vortex matter to defect 
lines in crystals. The procedure is applicable for sim- 
ulations as well as for experiments which provide real 
space structural information. We explained how a local 
cylindrical coordinate system can be defined on a coarse- 
grained scale along the line, in which the normal varia- 
tion with respect to the undulating line object can then 
be determined. The optimal length scale for the intro- 
duction of cylindrical reference frames turned out to be 
the persistence length of the fluctuating line object. We 
showed for the example of the density variation around 
the backbone of a bottle brush polymer that a normal 
variation around a flexible backbone determined for this 
length scale reproduces the normal variation around a 
rigid backbone on short and intermediate scales, where 
they should ideally be identical. This example is im- 
portant, since such stimuli-responsive polymers may find 
applications as sensors and actuators |30|, and for this 
purpose, as well as for asserting biomolecular functions 
of bottle-brushes such as aggrecane in mammalian joints 
[3l| , a deeper understanding of bottle-brush properties is 
required. 

For a reliable interpretation of experimental data on 
the variation of physical properties around an undulat- 
ing line object one has to go beyond the approximation 
of infinite rigidity. This manuscript suggests a controlled 
way to do this and to assess to what extent theoretical 
descriptions assuming infinite rigidity are valid. Besides 
the realm of synthetic and bio-polymers, we envisage po- 
tential applications of our concepts to problems such as 
the radial distribution of vacancies, solute atoms, etc., 
around lines in dislocation networks in crystals and liq- 
uid crystals. Also, heterogenous nucleation along line 
defects will be affected by their local curvature and it 
will be interesting to elucidate the differences to the case 
of completely rigid defect lines. 
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